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Proposed Work 


A research program was proposed for the testing and implementation of advanced turbulence 
models for non-equilibrium turbulent flows of aerodynamic importance that are of interest to 
NASA. Turbulence models that are being developed in connection with the Office of Naval 
Research ARI on Nonequilibrium Turbulence are provided for implementation and testing in 
aerodynamic flows at NASA Langley Research Center. Close interactions were established with 
researchers at NASA Langley RC and refinements to the models were made based on the results of 
these tests. The models that have been considered include two-equation models with an 
anisotropic eddy viscosity as well as full second-order closures. Three types of non-equilibrium 
corrections to the models have been considered in connection with the ARI on Nonequilibrium 
Turbulence : conducted for ONR 

(1) Anisotropies in the turbulent dissipation rate through an analysis of the transport equation for 
the tensor dissipation. The leading order contribution of this effect is through the addition of 
nonlinear strain dependent terms in the modeled scalar dissipation rate equation via the production 
coefficient C £l . The traditional constant value chosen for this coefficient makes it impossible to 

describe both equilibrium flows with moderate strain rates and non-equilibrium flows with large 
strain rates. 

(2) Non-equilibrium vortex stretching in the turbulent dissipation rate equation. The commonly 
used modeled transport equation for the turbulent dissipation rate is based on an equilibrium 
hypothesis whereby the production of dissipation by vortex stretching is exactly counter-balanced by 
the leading order part of the destruction of dissipation term. In order to describe departures from 
equilibrium, unbalanced vortex stretching will be allowed for which is described by a physically 
based relaxation model. 

(3) Non-equilibrium pressure-strain effects. Terms that are nonlinear in the mean velocity 
gradients are introduced into the model for the pressure-strain correlation through the 
implementation of a relaxation time approximation to a non-equilibrium algebraic stress model that 
bridges the equilibrium solution to the RDT solution for shear flows via a Padd approximation (see 
Appendix A). 

These models have the potential to lead to a new generation of Reynolds stress closures. 
While, as part of the Office of Naval Research ARI on Nonequilibrium Turbulence, the models 



will be tested in practical Naval Hydrodynamics flows, it would also be useful to test them in high 
speed aerodynamic flows that are of interest to NASA. This forms the raison d'etre of the present 
research. 

Research Accomplished 

The research focused on two central issues: 

(a) The development of a more robust regularization scheme for explicit algebraic stress 
models which form a cornerstone of the models being developed. The previously derived 
regularization scheme allowed the eddy viscosity to get too low when the mean strain rates became 
large. The new regularization procedure allows the eddy viscosity to approach a sufficiently large 
enough finite lower bound for numerical robustness. This has led to the better calculation of 
aerodynamic flows (see Appendix A). 

(b) The systematic incorporation of the effects of anisotropic dissipation into explicit 
algebraic stress models. By using the algebraic anisotropic dissipation rate model of Speziale and 
Gatski (1995), the explicit algebraic stress model approximation was repeated via integrity bases 
methods. It led to an explicit algebraic stress model where the coefficients simply assumed 
readjusted values (see Appendix B) . Thus the effects of anisotropic dissipation can be now be 
systematically implemented within the framework of a model that is only slightly more 
computationally expensive than the K - e Model. 

This research has great promise for future aerodynamic computations. Additional tests are 
currently underway. 
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Abstract 

An explicit algebraic stress equation, 
developed by Gatski and Speziale, is used in the 
framework of the K-e formulation to predict 
complex aerodynamic turbulent flows. The 
nonequilibrium effects are modeled through 
coefficients that depend nonlinearly on both 
rotational and irrotational strains. The 
proposed model was implemented in the ISAAC 
Navier-Stokes code. Comparisons with the 
experimental data are presented which clearly 
demonstrate that explicit algebraic stress 
models can predict the correct response to 
nonequilibrium flows. 

I. Introduction 

Computational fluid dynamics has become 
an increasingly powerful tool in the aerody- 
namic design of aerospace vehicles as a result of 
improvements in numerical algorithms and 
computer capabilities (e.g., speed, storage). 
Major future gains in efficiency are expected to 
come about as massively parallel supercomputer 
technology matures. However, some critical 
pacing items limit the effectiveness of computa- 
tional fluid dynamics in engineering. Chief 
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among these items is turbulence modeling. 
Numerous turbulence models of varying degrees 
of complexity, which can be classified as either 
eddy viscosity or full Reynolds stress models, 
have been proposed. Excellent reviews of turbu- 
lence models have been recently provided by 
both Speziale 1 and Wilcox. 2 

Eddy viscosity models use the Boussinesq 
isotropic effective viscosity concept, which as- 
sumes that the turbulent stresses in the mean 
momentum equation are equal to the product of 
an eddy viscosity and a mean strain rate. Zero-, 
one-, and two-equation models are among the 
most popular eddy viscosity models for 
engineering applications because of their ease of 
implementation in computational fluid dynam- 
ics codes. Algebraic or zero-equation models, 
which assume local equilibrium of the turbulent 
and mean flow, have provided reasonable 
predictions for simple flows. When the turbu- 
lent transport is important or the mean 
conditions change abruptly, these models do not 
work well. One-equation models improve the 
predictions for simple near-equilibrium flows 
but do not account for more complex effects on 
turbulence. Two-equation models are developed 
to take explicit account of the history of the 
turbulence through two transport equations for 
combinations of the turbulent length and time 
scales. These models offer good predictions of 
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the characteristics and physics of simple 
separated flows and flows with gradual changes 
in boundary conditions. However, basic two- 
equation models fail in many practical flows 
because they cannot properly account for 
streamline curvature, rotational strains and 
buoyancy; they provide an incorrect response to 
strong adverse pressure gradients; and they 
cannot describe the anisotropy of turbulence. As 
a result, various ad hoc modifications to these 
models have been proposed to achieve the 
proper response (see Lakshminarayana 3 ). In 
these modifications, effects on turbulence, such 
as those due to streamline curvature, have been 
directly accounted for in the eddy viscosity 
expression or have been reflected indirectly in 
the turbulence-model equations by modifying 
the dissipation-rate equation. The improved 
two-equation models predict a wider range of 
flows; however, they still fail to properly capture 
the physics in a broad class of flows. To 
overcome some of these deficiencies, two- 
equation turbulence models that are nonlinear 
in the mean strain rate were proposed by 
Speziale 4 and Rubinstein and Barton . 5 These 
models have provided accurate predictions of 
turbulence intensities. However, these models 
are not consistent with full Reynolds stress 
models because they have constant coefficients. 

Full Reynolds stress models represent the 
highest level of closure that is currently feasible 
for practical calculations. These models are 
superior to the two-equation models in that they 
eliminate the assumption that the turbulent 
stresses respond immediately to changes in the 
mean strain rate. Also, they account for the 
anisotropy of turbulence and body force effects 
on turbulence (e.g., due to streamline curvature 
and rotation) through extra production terms 
that explicitly appear in the Reynolds stress 
transport equation. However, models for many 
unknown turbulent quantities are required. 
This need is generally met by assuming that the 
turbulence is locally homogeneous and in equi- 
librium. Existing Reynolds stress models have 
been shown to give good descriptions of two- 
dimensional mean turbulent flows that are near 
equilibrium. However, computer costs and 
numerical stability problems that arise from the 


absence of a turbulent viscosity make assess- 
ments of the limitations of these models in 
predicting complex flows difficult. However, 
second-order closure models could be used to 
derive better two-equation models because 
fundamentally they are constructed on a 
stronger theoretical basis than the lower level 
models. 

Recently, a methodology for deriving a 
general nonlinear constitutive relation (or an 
explicit algebraic stress equation) for the 
Reynolds stress tensor from second-order 
closures, has been proposed by Gatski and 
Speziale , 6 based on the ideas of Pope . 7 This 
derivation is based on the assumptions that the 
net convection of the turbulent stresses is pro- 
portional to the net convection of the turbulent 
kinetic energy and that the structural parame- 
ters of the turbulence are constant along a 
streamline. As a result, a new generation of 
non-linear two-equation models is obtained with 
coefficients that depend on rotational and 
irrotational strains. This new feature extends 
the range of applicability of the standard two- 
equation models. 

Abid et al . 8 used the explicit algebraic 
stress relation within the context of the K-(o and 
K-e two-equation format to predict separated 
airfoil flows. The Launder, Reece and Rodi 9 
pressure-strain correlation model was consid- 
ered in the above study. Comparisons with the 
experimental data have shown that this new 
nonlinear turbulence model improves the ability 
of two-equation models to account for nonequi- 
librium effects. However, the Reynolds stress 
anisotropies were not well predicted. 

In this paper, the algebraic stress relation 
is applied within the context of the K-e two- 
equation format using the Speziale, Sarkar and 
Gatski 10 pressure-strain correlation model. The 
ability of the proposed model to predict complex 
flows which include nonequilibrium and 
anisotropic effects is assessed. Transonic flows 
over two airfoils and a wing are considered in 
this study. The ISAAC Navier-Stokes code is 
used to compute the three test cases. 
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II. Theoretical Analysis 

For a weakly compressible turbulent flow 
at high Rey nold s numbers, the Reynolds stress 
tensor = u t Uj is a solution of the transport 
equation 11 


Reynolds stress closures: the convection term 
minus the diffusion term in the Reynolds stress 
equation is proportional to the convection term 
minus the diffusion term in the turbulent 
kinetic energy equation and the Reynolds stress 
anisotropy 6^ is constant along a streamline. 
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The substitution of (3) and (4) into (1) 
yields the following algebraic stress equation: 
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given that riy is the pressure-strain correlation, 
D l is the turbulent transport term, e is the 
turbulent dissipation-rate, v is the kinematic 
viscosity, u) is the mean-velocity component, 
and p is the mean density. Explicit compress- 
ibility effects are neglected in Eq. (1) due to the 
applicability of Markovin’s hypothesis in these 
weakly compressible flows. 

If we contract the indices in (1), then we 
obtain the transport equation for the turbulent 
kinetic energy K = u,u, / 2 : 

= P-e + D% + vV 2 K (2) 

given that P = is the turbulence 
production term and Dj is the turbulent 
transport term. 

Rodi 12 proposed the idea of algebraic stress 
closure, which provides algebraic equations 
without solving differential equations for the 
Reynolds stresses. He assumed that 
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is the Reynolds stress anisotropy. Physically, 
two assumptions are made in the algebraic 
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where 


and 


are the mean-rate-of-strain tensor and mean- 
vorticity tensor, respectively. 

Given a pressure-strain-correlation model, 
(6) provides an implicit algebraic equation for 
the determination of the Reynolds stress x ir 
Computations that use this model have shown 
that stable numerical solutions can be difficult 
to obtain. Hence, an explicit algebraic stress 
equation which is a mathematically consistent 
representation of (6) is preferable. 

Pope 7 developed a methodology for obtain- 
ing explicit algebraic stress equations by using a 
tensorial polynomial expansion in the integrity 
basis. 5 Gatski and Speziale 6 used this method 
to derive an explicit algebraic stress equation for 
two- and three-dimensional turbulent flows. In 
order to generalize their results, they applied 
their algebraic stress representation to the 
general class of pressure-strain correlation 
models for which are linear in the 
anisotropic tensor 6y. The general linear form 
of H y is 
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The explicit nonlinear constitutive equation, 
derived by Gatski and Speziale, 6 is then given 
after regularization by 


which is equivalent to Eq. (12) to order r/ 4 and 
| 4 . Relation (17) does not change the value of 
near equilibrium conditions, but limits to 
a small non-zero value (= 0.2oq) for high values 
of 77 or £ to avoid numerical instabilities. In the 
present study, the pressure-strain-correlation 
model of Speziale, Sarkar, and Gatski 10 is 
considered; the coefficients are: 

C x = 6.8, C 2 = 0.36, C 3 = 1.25, 
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where p is the mean density and co = e/K is 
the specific dissipation rate. The constants in 
(11W13) are given by 

«x=(f-C 2 )f, a 2 =(2-C 3 ) 2 ^ (14) 

«3=(2-C 4 ) 2 ^-, (15) 


a 5 =(2-C 3 )g, g = 
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(16) 


To avoid numerical problems in the initial 
stages of the computation or in the free-stream 
region, a modified form of C* is used 
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C 4 = 0.40, C 8 = 1.88 


(18) 


The nonlinear constitutive equation (10) 
must be solved in conjunction with the following 
modeled transport equations. 
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given that p t( =C' fll p and C^(= 0.081) is the 

value of Cl in the fogarithmic layer. The 
coefficients of the model are 

<j k = 1.0, k = 0.40, C t2 = 1.83, C Ci = 1.44 
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given that u T is the shear velocity and y is 
normal to the wall. Note that new model can be 
integrated directly to the wall without adding a 
damping to the eddy viscosity. The function f is 
introduced to remove the singularity in the 
dissipation rate equation at the wall. 

At the wall, the boundary conditions for K 
and e are 
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m. Results and Discussion 

The calculations to be presented were done 
with the three-dimensional Navier-Stokes 
ISAAC code, 13 which uses a second-order accu- 
rate finite-volume scheme. The convective 
terms are discretized with an upwind scheme 
that is based on Roe's flux-difference splitting 
method. All viscous terms are centrally differ- 
enced. The equations are integrated in time 
with an implicit, spatially split approximate- 
factorization scheme. 

The performance of the explicit algebraic 
turbulence model (hereafter referred to as 
EASM) was evaluated for the flat-plate turbu- 
lent boundary layer at a zero-pressure gradient. 
As expected (the results are not shown here), 
the turbulence model yielded good predictions 
for the mean-velocity profiles and skin-friction 
coefficients. Although some turbulence proper- 
ties near the wall are not captured (i.e., the peak 
of the turbulent kinetic energy), the algebraic 
stress model does give accurate results away 
from the buffer layer (i.e., y + > 30). Remember 
that the algebraic stress model can be inte- 
grated directly to a solid boundary with no 
damping function in the turbulent eddy 
viscosity. 

The first two test cases to be considered are 
the RAE 2822 airfoil flows (cases 9 and 10), 
which were tested by Cooke et al. 14 The airfoil 
has a maximum thickness of 12.1 percent c and 
a leading-edge radius of 0.827 percent c (c is the 
chord of the airfoil). The grid used is a 257x97 C 
mesh with 177 points on the airfoil, and a 
minimum spacing at the wall of 0.932xl0^c. 
The outer boundary extent is approximately 18c, 
and transition is assumed at 3 percent c. For 
the case 9, the conditions include a Mach 
number M <*, = 0.73, an angle of attack a = 2.8°, 
and a Reynolds number Re = 6.5xl0 6 . This case 
contains no separated flow. For the case 10, the 
conditions include a Mach number Moo = 0.75, 
an angle of attack a = 2.72, and a Reynolds 


number Re = 6.2xl0 6 . This case involves sepa- 
ration based on visual surface streamline 
patterns. However, there are no skin-friction 
coefficient data indicating separation. Hence, 
case 10 is considered as an incipiently separated 
flow and, therefore, is more challenging than the 
previous case. 

Figures 1 and 2 compare the surface 
pressure and skin-friction coefficients computed 
along the airfoil surface with the experimental 
data for case 9. It is clear that the explicit alge- 
braic stress model provides a good representa- 
tion of the pressure over most of the airfoil. 
However, the turbulence model over predicts the 
skin-friction coefficient downstream of the 
shock. This deficiency results from the tendency 
of the models based on K-e formulation to 
predict excessive near-wall levels of turbulent 
length scale in the presence of an adverse pres- 
sure gradient, which leads to high values of the 
eddy-viscosity. A modification of the dissipation 
equation is required in order to improve the 
response of the algebraic stress model to adverse 
pressure-gradient effects. 

In order to demonstrate the improvement 
resulting from the use of the EASM model for 
non-equilibrium flows, comparisons between the 
results obtained by the EASM model and the 
Speziale, Abid and Anderson K-e model 15 
(hereafter referred to as SAA) were performed 
(Figures 3-10). From Figure 4, it appears 
clearly that neither turbulence model predicts 
separation. This is reflected by the high level of 
the skin-friction coefficient, downstream from 
the shock. This probably is a result of the 
inability of the length scale equation to provide 
proper response to adverse pressure gradients. 
To date, several modifications to the dissipation 
equation for separation do not seem to be 
successful. On the other hand, the EASM model 
predicts the shock location better than the SAA 
model, although slightly downstream of the 
experimental shock location (see Figure 3). This 
results from the prediction by the EASM of 
lower values of eddy viscosity in the inner part 
of the boundary layer, therefore, lower values of 
the turbulent kinetic energy (see Figure 7). 
Comparison of the computed and measured 
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velocity profiles further support the latter 
observation. An additional finding that can be 
inferred from the above comparison is that the 
EASM model gives a realistic representation of 
the normal stresses (see Figures 8-10). 

The third test case to be considered is the 
ONERA M6 wing at Mach number of 0.8447, an 
angle of attack a of 5.06 and a Reynolds number 
of 11.7xl0 6 based on the mean aerodynamic 
chord. 16 A C-0 grid, used in this study has 
193x49x33 points in the streamwise, normal 
and spanwise direction. The minimum normal 
spacing over the wing of 0.000015 c roo t and a 
distance from the wing to the outer boundary of 
at least 7.95 c roo t . No wind tunnel test correc- 
tions are employed for this case. 

Figure 11 shows a comparison of the 
surface pressure distributions with the experi- 
mental data at four different spanwise locations 
2 y!B. It is clear from this figure, that the 
predicted shock location and the surface 
pressure distributions by the EASM model are 
in good agreement with the experimental data, 
and similar to the results reported in [17] for the 
Johnson-King model, which has been highly 
tuned for airfoil flows. 

Conclusions 

A study of an explicit algebraic stress 
model, used in the framework of the K-e 
formulation for separated turbulent flows, has 
been conducted. This new generation of two- 
equation models, which is derived from second- 
order closures, has been tested against three 
test cases, two of which involve separation. Two 
major findings have been made in this study: 
explicit algebraic stress models have shown 
some improvement over the standard two- 
equation models because of their ability to 
account for nonequilibrium effects and to give a 
realistic representation of the anisotropy of the 
turbulence. However, this improvement is still 
limited by the dissipation rate equation which 
fails to respond properly to adverse pressure 
gradients. A major research effort to correct 
this deficiency is currently underway. 
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Figure 7: Comparison of turbulent kinetic energy distributions for RAE 2822 airfoil 
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Figure 8: Comparison of u 2 normal stress distributions for RAE 2822 airfoil (Case 10) 
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Figure 9: Comparison of v 2 normal stress distributions for RAE 2822 airfoil (Case 10) 



Figure 10: Comparison of w 2 normal stress distributions for RAE 2822 airfoil (Case 10) 
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1. INTRODUCTION 

Turbulent flows near solid boundaries - or at low turbulence Reynolds numbers - can 
exhibit significant anisotropies in the turbulent dissipation rate 1 . Nevertheless, Reynolds 
stress turbulence closures are routinely formulated that neglect such effects by invoking the 
Kolmogorov assumption of local isotropy 2 . Recently, however, attempts have been made 
to extend full Reynolds stress turbulence closures to incorporate the effects of anisotropic 
dissipation (see Speziale, Raj and Gatski 3 , Speziale and Gatski 4 , Oberlack 5 and Hallback et 
a/. 6 ). These more sophisticated Reynolds stress turbulence closures can involve the solution 
of up to twelve additional transport equations. As such, most of these models are not 
currently feasible for the solution of complex turbulent flows in an engineering setting. 

During the past few years, explicit algebraic stress models have been developed that are 
formally consistent with full second-order closures in the limit of homogeneous turbulence in 
equilibrium (see Gatski and Speziale 7 ). These models allow for the solution of complex tur- 
bulent flows with a substantially reduced level of computation compared to full second-order 
closures, since they constitute two-equation models 7,8 . The purpose of the present note is to 
show how the effects of anisotropic dissipation can be systematically incorporated into these 
explicit algebraic stress models by a simple readjustment of the coefficients. For homoge- 
neous turbulent flows that are close to equilibrium, it will be shown that the results obtained 
from such models are virtually indistinguishable from those obtained from a full second-order 
closure model with the anisotropic dissipation rate model of Speziale and Gatski 4 . All of 
this extra turbulence physics is incorporated within the framework of a two-equation model 
that is not much more computationally expensive to implement than the standard K — e 
model. 
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2. THEORETICAL BACKGROUND 

We will consider incompressible turbulent flows where the velocity u, and kinematic 
pressure P are decomposed, respectively, into ensemble mean and fluctuating parts as follows: 

Vi = Vi +U{, P - P + p. ( 1 ) 

In homogeneous turbulence, where all higher-order correlations are spatially uniform, the 
Reynolds stress tensor r t j = v~Uj satisfies the transport equation 9 

dvj dvi 

T ij = r lk — — Tj k h n,j — £ij , (2) 


where 

n _ / du i , du J x 

_ o ® U ' ® u i 
£lJ ~ V dx k dx k 

are, respectively, the pressure-strain correlation and the dissipation rate tensor. Thus, in 
homogeneous turbulence, only n,j and £,j need to be modeled in order to achieve closure. 

Speziale, Sarkar and Gatski 10 showed that, for two-dimensional mean turbulent flows in 
equilibrium, the commonly used hierarchy of pressure-strain models simplifies to: 

II,j = —C\ebij + C 2 e(bi k b k j — ^b k (b k ( 8 ij) 

-\-C 3 KSij + C 4 K(b lk Sj k + bj k Si k (3) 


where 


3 b k iS k i 8 ij) T C b K{b ik u jk T bj k uj { k ) 


q.. = h^ 4 -^ n-. - h - Til Z 3 j 

13 2 ^dxj + dx/’ ,3 2^dxj dxi ,J 2 K 


I KSij 


are, respectively, the mean rate of strain tensor, the mean vorticity tensor, and the Reynolds 
stress anisotropy tensor; K = is the turbulent kinetic energy. The SSG model is a 
simple extension of (3) that is valid for moderate departures from equilibrium. It has been 
found that the nonlinear return term containing C 2 can be neglected without introducing an 
appreciable error. 7,10 With the choice of constants 


C x = 6.80, C 2 = 0, C 3 = 0.36, C 4 = 1.25, C 5 = 0.40, 


in (3), excellent equilibrium values are obtained for the benchmark case of homogeneous 
shear flow. We refer to this as the linearized, equilibrium form of the SSG model. 
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The Kolmogorov assumption of local isotropy is typically invoked wherein it is assumed 
that 2,9 

2 

e *i = g (4) 


where e = |e„ is the (scalar) turbulent dissipation rate. It is generally accepted that 
homogeneous turbulent flows, with constant mean velocity gradients, achieve equilibrium 
values for b tJ that are largely independent of the initial conditions. This is characterized by 




= 0 


or, equivalently, 

r,i = 2(v - e)i, :y + \(T - e)6 h (5) 

where V = —Tijdvi/dxj is the turbulence production. The substitution of (3) - (5) into 
(2) yields an implicit algebraic system that can be solved by integrity bases methods. This 
solution - which has come to be referred to as an explicit algebraic stress model (ASM) - is 
given in the equilibrium limit of homogeneous turbulence by: 7 




2„ t 3 f 

3 S,i PE 5 " 


where 


3 

+<* 1 — ; + SjkVki) 


-a 2 


K 3 (- _ 1- _ „ 

ySikSkj — -S kl S kl 8, 


•)] 


a o — (jj — C 3 ) 9 , ai — 2 (3 _ ^ 3 ) ^ ~ C*,)g‘ 
a 2 = (g — C 3 ) (2 — Ci)g 2 , g = (-Ci + — — l) 
v = —7 ( 3«$i) 1/2 , i = ~ fei) 1/2 

l ao £ OLQ E 


-1 


( 6 ) 


( 7 ) 


and V/e assumes its constant equilibrium value. When far from equilibrium, a singularity 
may arise since the denominator (3 — 2r] 2 + 6£ 2 ) in (6) can vanish for sufficiently high strain 
rates 77. Gatski and Speziale 7 introduced a regularized expression for 3/(3 — 27y 2 + 6f) 
which eliminated the singular behavior. However, that model is not formally valid for non- 
equilibrium turbulence — particularly in the rapid distortion limit. 
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Speziale and Xu 11 later introduced a formal Pade approximation that built in some 
limited consistency with Rapid Distortion Theory (RDT) for homogeneous shear flow. They 
rewrote (6) in the form: 

1^ ^ 


2 2 3 / j V 

r \j — — QJq - I ~ “1“ ° i 2~ 2~ ySikSfcj ^ $ kl Skl^ij J (8) 


~ *J ”1 2 IK—KJ I KIJ l — 2 2 

c c c 


and made use of the fact that in the short-time RDT solution (77 — > 00), Kj K 0 remains of 
order one. This implies that 

<=i ~ i. (9) 

V 

It is obvious that the equilibrium model (6) violates this constraint (ag ~ l/ 7 ? 2 instead). 
Speziale and Xu 11 then introduced a Pade approximation and obtained the expression: 

_ (1 + 2f 2 )(l + 6 t? 5 ) + 

0 (l + 2e)(l + 2£ 2 +y + 6A V «) 0 1 1 

(with the constant / 3 0 7) to ensure asymptotic consistency and the proper energy growth 

rate in line with the RDT data. By a comparable Pade approximation, they also derived 
the expression 

P+jgKi+jVjsL .. (11) 

' (i+2p)(i+2e+ftii 6 ) 1 ’ 

for i — 1,2 (with « 6 and /? 2 ~ 4), in order to establish consistency with the approach to a 
one component turbulence in the RDT limit of homogeneous shear flow. For near-equilibrium 
turbulent flows, (8) with (10) - (11) yields results that are virtually indistinguishable from 
( 6 ). 

3. ANISOTROPIC DISSIPATION 

In a recent study, Speziale and Gatski 4 derived a modeled transport equation for €ij 
which is valid for homogeneous turbulence. By invoking the equilibrium limit where 


dij = 0 


( 12 ) 


for the anisotropy of dissipation they obtained an algebraic system — analogous to that 
for algebraic stress models — which was solved by integrity bases methods. This ultimately 
led to the algebraic model: 4 

iij = -2C„,{rS„ - + S jk m ki ) 


( 30 n _ 2 W 2 
Ml® 3 11 ) T 7? 


1^ ^ 


+ c e lfv/e- l {SikSkj ~ 3 s ~n s ~«Sii)} 


(13) 
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where 


c • • — —eS • • 
d = £o 3^ 

J 2 e 


(14) 


is the anisotropy of dissipation and r = K/e is the turbulent time scale. Here, 

1 


G{it 15(^5 +V/e-l ) {1 + 2r ‘ 


_ - T 2 ( n a3 i i 
3 \ C e 5 + V/e 


( rr a 3 + rr V 

^s + P/e-lj 

y 






(15) 


and 


C’es — 5.8, <23 — 0.6 

are constants. This then yields the full dissipation rate tensor since, by definition, 

2 


'•j 


e6^,j “H 2edij. 


(16) 


After introducing the anisotropic dissipation model for in (2), we can derive an alge- 
braic stress model with anisotropic dissipation. Again, making use of the fact that b t j = 0 
for equilibrium turbulent flows, the Reynolds stress transport equation (2) then reduces to: 

2(P — e)bij = —2K{b lk Sjk + b ]k S lk — - b mn S mn 8 ij ) — - KSij 

—2K(bikuJjk + bjkUJik ) + n tJ — 2 edij. (17) 

The explicit ASM incorporating anisotropic dissipation - which is obtained from the solution 
of (17) after (3) and (13) are implemented - is of the same general tensorial form as (6): 

b,j = G {1) T[p + G {2] Tp ] + G {3) T[P 






(18) 


where 


T-p = Sij , T-p = S ik u kj + S jk u ki 

Tip = S tk S kj - ^S mn S mn 6 tJ 


are the integrity bases. The solution is given by 


<3 (1) = -ao T 


~ 1 — Ai + \tj 2 A 3 + 2fA 2 


1 - h 2 + 2e 2 


G ( 2 ) 


-air 


— 1 ~ Ai — A 2 + 3^/ 2 ( A3 + 2A 2 ) 


1 - w + W 
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G< 3 > = ia 2 r J 
4 


2-A 3 + 2A 1 - 2 e{A 3 + 2A 2 ) 


1 - iv 2 + 2£ 2 


( 19 ) 


where 


A? = 


— 


Ar = 

2 


4a 


fJC 


i-c 3 

4fi*(n a 3 + n) 


s(l - C 3 )(2 - C s ) C.s + 2-1 

2 4C„,(f?a 3 - £) 


ff(| - Cj)(2 - C 4 ) CU + ?-l 


Equivalently, from (18) we have 


Tij = \K8 tj + 2K{G [x) Tl- ] + + G {3) t[- ] ) 


( 20 ) 


which is obviously of the same tensorial form as (6) — only the coefficients are different. 

The factor 1/(1 — | tj 2 -f 2£ 2 ) in (19) can be regularized in the same general way as dis- 
cussed earlier to ensure the correct asymptotic behavior in the RDT limit. The standard 
explicit ASM given in (6) is then recovered in the limit as C ^ (and hence, A 2 , and A 3 ) 
— > 0. Anisotropies in the dissipation rate are then accounted for simply through a systematic 
readjustment of the coefficients. 

4. RESULTS AND DISCUSSION 

The anisotropic dissipation rate model has been tested in detail by Speziale and Gatski 4 
within the context of a full second-order closure, which will not be repeated here. Our 
purpose in this note is to simply demonstrate that - for homogeneous turbulent flows close 
to equilibrium - the new explicit ASM with anisotropic dissipation yields results that are 
indistinguishable from the full second-order closure, with anisotropic dissipation, on which 
it is based. 

The new explicit ASM with anisotropic dissipation derived herein is solved with a modeled 
transport equation for the turbulent dissipation rate that is of the form 4 

i = c; x ^-v-cA (2i) 


for homogeneous turbulence. Here 

2(1 + a) 


c ; ! = 1.26 + 


C t5 + 2C li rj 2 - 1 


15 [{C th + 2C #J r f - l) 2 + 2/? 2 £ 2 - |/? 2 VJ 


( 22 ) 


3,14 16. 

a = — ( — a-i ) 

j 00 / 


33 


where 
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a _ 7 , 1 * 15 1 

A “TI Q 3 + n’ = 


and 0^ and C t2 are constants that assume the approximate values of 0.09 and 1.83, respec- 
tively. In (22), T) = (SijSijY^T and £ = The effects of anisotropic dissipation 

are rigorously accounted for in (21) through the variable coefficient C* x . For inhomogeneous 
turbulent flows, a gradient transport term of the standard form 


vj de . 

<r t dxi 

is added to the r.h.s. of (21) to account for turbulent diffusion. 

In order to demonstrate the ability of the new model to properly capture the physics of 
the more complicated full second-order closure with anisotropic dissipation, we will present 
a simple example of a benchmark turbulent flow. In Table 1, we compare the equilibrium 
results of the new explicit ASM for homogeneous shear flow with the predictions of the 
linearized SSG second-order closure incorporating the anisotropic dissipation rate model of 
Speziale and Gatski 4 . These calculations were conducted with a fourth-order accurate Runge- 
Kutta numerical integration scheme. It is clear from these calculations that the new explicit 
ASM yields results that are virtually identical to the results of the full second-order closure 
with anisotropic dissipation and compare favorably with the DNS results of Rogers, Moin 
and Reynolds 12 . This definitively demonstrates the power of the new model to yield results 
that are indistinguishable from a full second-order closure — with anisotropic dissipation 
— for turbulent flows that are close to equilibrium. All of this at a small fraction of the 
computer costs. While the effects of anisotropic dissipation are not that significant for this 
case, they can be important in other flows of engineering interest as alluded to earlier. We 
consider this to be a highly promising new approach for such flows that can have important 
engineering applications and, thus, warrants further study. 
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Equilibrium 

Values 

Explicit 

ASM 

Full 

Closure 

DNS 

b n 

0.204 

0.205 

0.215 

b\2 

-0.150 

-0.150 

-0.158 

b 22 

-0.148 

-0.147 

-0.153 

&33 

-0.056 

-0.058 

-0.062 

SK/e 

5.98 

5.96 

5.70 


Table 1. Equilibrium values for homogeneous shear flow. Comparison of the new explicit 
ASM incorporating anisotropic dissipation with the DNS of Rogers, Moin and Reynolds 12 
and a full Reynolds stress closure with anisotropic dissipation (containing the linearized SSG 
second-order closure and the anisotropic dissipation rate model of Speziale and Gatski 4 ). 




